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Abstract. We compare different partitioning schemes for the box-counting algorithm in the multifractal 
analysis by computing the singularity spectrum and the distribution of the box probabilities. As model 
system we use the Anderson model of localization in two and three dimensions. We show that a partitioning 
scheme which includes unrestricted values of the box size and an average over all box origins leads to smaller 
error bounds than the standard method using only integer ratios of the linear system size and the box size 
which was found by Rodriguez et al. [T] to yield the most reliable results. 



considered in the comparative study by Rodriguez et al. 
Therefore, we compare in this paper the fit quality of the 
best method proposed in [I] with the participation scheme 
allowing non-integer ratios N /L [3] . We show that the lat- 
ter partitioning scheme yields smaller errors in the least 
squares fits than the standard scheme which allows only 
integer ratios of N /L. 

The paper is structured in the following way; First, 
we briefiy introduce the Anderson model in Sec. [2] and 
the multifractal analysis in Sec. [31 We then discuss differ- 
ent partitioning schemes for the box-counting algorithm 
in Sec. 0] and compare their performances in Sec. [5] This 
is followed by a brief conclusion in Sec. [51 



1 Introduction 

Systems with multifractal properties are often found in 
nature, e.g. in biological systems, in meteorology and ge- 
ology, or in economics, e.g. stock market time series. For 
a detailed study of these properties multifractal analysis 
plays an important role. Also the electronic wave func- 
tions 'P of physical models like the Anderson model of lo- 
calization at the mobility edge, quasiperiodic tilings and 
chaotic systems show multifractal properties [2131415] . One 
can distinguish two methods for the multifractal analysis: 
the system-size scaling approach and the box-size scaling 
approach |^ . In both approaches the system is partitioned 
into boxes and one measures the probability to find the 
electron in a particular box, i.e., we study the spatial dis- 
tribution of the wave functions T . While for the system- 
size scaling the box size is kept constant and the system 
size is varied, for box-size scaling one studies how the prob- 
ability distribution changes with the box size. The latter 
method is more common because it requires less comput- 
ing time and storage. 

For box-size scaling different partitioning schemes are 
known, where the scheme has a crucial influence on the 
fit results and the statistical analysis. Rodriguez et al. 
compared several approaches for partitioning the system 
into boxes including different box shapes and box sizes as 
well as an averaging over different box origins |1I6] . They 
found that a partitioning of the system into cubic boxes 
with integer ratios of the linear system size N and the lin- 
ear box size L is numerically most reliable in comparison 
to the other studied methods. However, a previous work 
by one of us showed that the fit quality can be improved 
by adapting the multifractal analysis to work also for non- 
integer ratios N / L and by averaging over different origins 
of the boxes [31819] . Unfortunately, this method was not 



2 Anderson Model of Localization 

As model system for studying the influence of the parti- 
tioning in the multifractal analysis we use the eigenstates 
of the tight-binding Anderson Hamiltonian 
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represented in the orthogonal basis of states |r) associated 
to a vertex r. The onsite potential is chosen accord- 
ing to a uniform distribution in the interval \—w/2,w/2\ 
for a given disorder strength w. Each site r is connected 
to the nearest neighbors in a regular square/cubic lat- 
tice with periodic boundary conditions. For a given re- 
alization of disorder we consider only the eigenstate #r 
of the Anderson Hamiltonian closest to the band center 
at -E = 0. The wave functions are obtained using the pro- 
gram JADAMILU, which is based on the Jacobi-Davidson 
method [TU] and allows one to compute selected eigenval- 
ues and eigenfunctions of large sparse matrices with high 
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Fig. 1: Probability amplitude |<?r| of a wave function on 
a square lattice of iV^ = 360 x 360 sites. Every site 
r with an amplitude larger than the average is shown, 
i.e., |«?r| > 1/-^ and four different grey levels distinguish 
whether |<l>r|^ > 2^ /N^ for j = 0,1,2,3. 



efficiency. In Fig. [T] we visualize a two-dimensional wave 
function for the disorder strength u; = 4 for the linear 
system size N = 360. 

The investigations in this paper are restricted to mod- 
erate system sizes with at most N — 360 in two dimen- 
sions and TV = 60 in three dimensions which are suffi- 
ciently large for the presented purpose. We have also stud- 
ied larger system sizes up to iV = 720 in two dimensions 
and up to TV = 90 in three dimensions but we have not 
found any significant differences for the fit results for the 
single wave functions. Of course, the system size has a cru- 
cial influence on the ensemble average, and the below dis- 
cussed measures are not yet completely size-independent 
for these system sizes. Therefore, large system sizes should 
be preferred. As we are only interested in the optimization 
of the partitioning scheme we mainly consider single wave 
functions. 



3 Multifractal Analysis 

The multifractal analysis is based on a standard box- 
counting algorithm, i.e., in d dimensions the system of 
size V = N"^ is divided into B — N'^/L'^ boxes of linear 
size L [2I5I11I12] . For each box b we measure the proba- 
bility to find the electron in this particular box. Hence, 
for an eigenstate <l>r and the box size L we determine the 
probability for the 6th box 5 



and construct for each parameter q a measure by normal- 
izing the qth moments of this probability 



P{q,<l>,L) 



where 



P(g,<P,i)=^M^(<Z',i) 



6=1 



The mass exponent of a wave function 
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can be obtained by a least squares fit of InP versus Ine, 
where e = L/N. The generalized dimensions are easily 
obtained as Dq — Tq/{q— 1). 

According to Eq. ([5]) small boxes should be considered 
for the least squares fit. Further, the power-law behavior 
defined in Eq. ([5]) requires the absence of a typical length 
scale in the distribution of the probability density of the 
wave functions. Therefore, it should be computed for box 
sizes I L N to properly display the multifractal prop- 
erties, where I corresponds to microscopic length scales of 
the system like e.g. the lattice spacing [T3]. Rodriguez et 
al. used box sizes from L = 10 to A^/2 [1]. For reasons 
of comparison we have performed most of the calculations 
also for this range, but we also show that the range of con- 
sidered box sizes has a crucial influence on the observed 
results. 

Another quantity often used for the multifractal anal- 
ysis is the singularity spectrum /(a^), which denotes the 
fractal dimension of the set of points where the wave- 
function intensity behaves according to \<Pr\ OC TV-"", 
which means that in a discrete system the number of 
such points scales as N^^°'i'> [T3]. The singularity spec- 
trum f{aq) follows from a parametric representation of 
the singularity strength 



lim 

£->0 
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with 



A{q,$,L) = ^^b(q,<?,L)ln/X;,(l,<?>,L) 
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and the fractal dimension 
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rG box b 



Both quantities can be determined by least squares fits 
of A{q, <!>,£) or F{q,(I>,L) versus Ine considering differ- 
ent box sizes L. The mass exponents Tq and the gener- 
alized dimensions Dq can also be obtained from the sin- 
(2^ gularity spectrum f{aq) via the Legendre transformation 
Tq = Dq{q - 1) = qaq - f{aq) [H]- 
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Regarding the partitioning of the system into boxes 
several options are possible, which are discussed in more 
detail in the following two sections. 



4 Partitioning Schemes 

In this section we briefly introduce the considered meth- 
ods for the optimization of the partitioning in multifractal 
analysis. In particular, we discuss the following three par- 
titioning schemes for the box-counting algorithm: 

A) cubic boxes with integer ratios N /L, 

B) cubic boxes with unrestricted values for L, and 

C) cubic boxes with integer ratios N / L including an av- 
erage over all box origins. 

To be specific, we define the left upper corners of a box as 
box origin in the figures below. 



4.1 Scheme A: Cubic boxes with integer ratios N/L 

This quite simple partitioning scheme is the most straight- 
forward procedure in which only such box sizes are consid- 
ered for which the boxes cover the system gapless without 
overlap as visualized in Fig. [51 Hence, the ratio 1/e = N/L 
is restricted to integer values. Obviously this condition is 
a severe limitation of the possible choices for N and L. For 
example, in the case of TV = 6 the two partitionings with 
L = 2 and L — 3 shown in Fig.[2]are the only possibilities. 
For TV = 5 or TV = 7 there does not exist any such way of 
subdividing the system. 

Exploiting periodic boundary conditions, one could uti- 
lize different box origins. But usually no additional aver- 
age over different box origins is considered. Rodriguez et 
al. studied several partitioning schemes and found that 
this one shows the most reliable performance of the con- 
sidered schemes T. However, we show in this paper that a 
method considering unrestricted values of box sizes shows 
a better performance. This scheme was already success- 
fully applied for the multifractal analysis I3J8] . Unfortu- 
nately, it was not taken into account in the comparison by 
Rodriguez et al. [I]. 





4.2 Scheme B: Cubic boxes with unrestricted values 
for L 

The fitting of the scaling exponents in the box-size scal- 
ing can be enhanced if all possible box sizes L as well 
as an average over all box origins is considered |3I8I9| . 
This method has a rather simple geometrical interpreta- 
tion, which is displayed in Fig. [31 First, the original sys- 
tem is periodically repeated LCM(L,A^)/A^ times (LCM 
labels the least common multiple) in each direction. The 
new system can then be covered gapless for any linear box 
size L. If N and L are coprime, then this method automat- 
ically includes an average over all different box origins as 
shown in Fig. [31 For the other cases we explicitly consider 
an average over all different box origins. 

From the computational perspective it is easier not 
to look at the periodically repeated system but rather to 
imagine that all boxes with origins outside the original 
system are folded back so that their origins fall into the 
original system. One can easily see from Fig. [3] that if TV 
and L are not coprime such a folding back would not yield 
all sites of the original system, but it is of course straight- 
forward to generalize the idea and to place a box of linear 
size L on each site of the original system in this case, too. 




Fig. 2: Partitioning scheme with integer ratio N/L for a 
system of linear size TV = 6 and box sizes L ~ 2 (left) and 
L — 3 (right). The red dots indicate the considered box 
origins. 



LCM(L,TV) 



Fig. 3: Partitioning scheme for unrestricted values of 
the box size L for a system of linear size TV = 4 and 
the box size L — 3. The original system is repeated 
LCM(i, TV)/TV — 3 times in each direction so that the re- 
sulting system of linear size LCM(i, TV) = 12 can be filled 
gapless by the boxes. For coprime L and TV this method 
intrinsically includes an averaging over all box origins as 
indicated by the red box origins. The arrows visualize this 
for the sites in the third column of the original system. 
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Hence, in all cases each of the N"^ sites is considered 
times and the proper average for a quantity X is obtained 

by 

W-;^ E ^- (10) 

box origins 

This average then replaces the simple sum over all B 
boxes. 

Rodriguez et al. also studied a method which allows 
one to utilize all box sizes L [T] . In contrast to our method, 
they did not periodically repeat the system LCM(L, N)/N 
times. Hence, they faced the problem that different effec- 
tive system sizes occur, and they had to introduce an ad- 
ditional average over these different effective system sizes. 
Our method does not require such an additional averaging 
because the boxes fill the new system gapless, so that no 
effective system sizes appear. 

4.3 Scheme C: Cubic boxes with integer ratios N/L 
including an average over all box origins 

This method is a mixture of the two former partitioning 
schemes. Now we only consider integer ratios but include 
the average over all box origins as in the previous method. 
Rodriguez et al. also studied the average over all initial 
box origins but pointed out that for an ensemble aver- 
age over many wave functions this method reduces the 
standard deviation almost to zero [T] . The reason for this 
apparently negative side effect is that they counted each 
box origin as a different disorder configuration, which is 
clearly not the case. 

Hence, a proper ensemble average over different disor- 
der configurations can still be obtained by computing the 
average quantities P(g, ^, L), A{q, <P, L) and F{q, L) for 
all box origins for each wave function according to Eq. (jlOD 
but without considering each box origin as a new disorder 
configuration. In this case no unwanted reduction of the 
standard deviation of the ensemble average occurs. How- 
ever, the average over the box origins in Eq. (jlOp means 
that the data for each wave functions are evaluated with a 
better accuracy, so that a reasonable reduction of the stan- 
dard deviation of the ensemble average can be expected. 

To avoid a misunderstanding, we point out that the 
replicas in Fig. |3] are used for an easy explanation of our 
scheme only. They are not necessary for our computation. 
One can also utilize the following interpretation: Due to 
periodic boundary conditions, no site of the system can or 
should be distinguished or preferred over any other site. 
But the choice of the system origin makes such a distinc- 
tion. This bias can be avoided by treating every other site 
equivalently, by considering every site as a possible sys- 
tem origin. Then the generic placement of boxes starting 
at all the possible system origins is equivalent to our above 
placement of box origins on each site of the system. 

5 Comparison of Results 

For a comparison of the three partitioning schemes we 
study the singularity spectrum of wave functions closest to 



the band center E ^ Q for certain disorder configurations. 
We consider wave functions in different dimensions and 
study also the influence of the considered range of box 
sizes. 



5.1 Two-dimensional wave functions 

Figures |4] and [5] show the scaling behavior of the quantities 
A{q,(I>,L) and F{q,<P,L) with the box size L for a two- 
dimensional wave function for the system sizes N — 350 
and iV = 360. The results can be well described by a 
power law. The stronger fluctuations for negative values 
of q are well known and can be explained by the less ac- 
curate values of the wave functions at vertices with small 
amplitudes, which are enhanced for negative powers of fit 
in Eq. Apart from fluctuations for very small box sizes 
and negative q, the data are sufficiently close to straight 
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Fig. 4: Scaling behavior of the quantities A{q, <P, L) (top) 
and F{q, <P, L) (bottom) with the ratio e = L/N of the box 
size L and the linear system size TV for the two-dimensional 
wave function of the Anderson Hamiltonian shown in Fig. 
[T] with N = 360 and disorder strength w — 4. The re- 
sults are shown for the scheme A with integer ratios N/ L 
(squares) and scheme B with unrestricted values of L and 
a box-origin average (circles) for 2 < L < N/2. The dot- 
ted lines correspond to least squares fits considering all 
box sizes in the range 5 < L < N/10 marked by the two 
vertical dashed lines. 
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Fig. 5: Same as Fig. 01 but for the system size N = 350. 

lines as visualized in Figs.|3]and[5]so that we can determine 
the singularity strengths aq {<P) and the fractal dimensions 
/g(^) from the slopes. If not stated differently we consider 
box sizes in the range from L = 10 to N/2 for the com- 
putation of the singularity spectrum in accordance to the 
setup used by Rodriguez et al. [1114) . 

The corresponding singularity spectrum f{aq) is shown 
in Fig. ini From literature it is known that f{aq) is a con- 
vex function of a^, where the maximum of f{oLq) corre- 
sponds to g = with /(ao) = d being the support of the 
measure |llll2j . This is verified in Fig. [S] with /(ao) ~ 2. 
Further, for q = 1 we have /(ai) = oti and the limits 
q — > ±oo yield the minimal/ maximal value of a^, which 
corresponds to the singularity associated to the box con- 
taining the largest/smallest measure |5]. For the periodic 
system, i.e. for vanishing disorder, the wave functions are 
Bloch waves extending over the entire system and there- 
fore the singularity spectrum is expected to converge to 
the point f{d) — d. In the limit of a strongly disordered 
system with localized wave functions, the spectrum con- 
verges to the points /(O) — and /(oo) = d. The wave 
functions of the two-dimensional Anderson Hamiltonian 
show large spatial fluctuations, which increase with in- 
creasing disorder strength 15! . However, for moderate dis- 
order strengths the wave functions can be approximately 
described by the multifractal analysis. 

Comparing the results for the three different partition- 
ing schemes we observe that the least squares fit of the 
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Fig. 6: Singularity spectrum f{aq) for the three different 
methods for the systems shown in Fig. 2] with N = 360 
(top) and in Fig. [5] with N — 350 (bottom) with consid- 
ered box sizes from L = 10 to N/2. The plot shows Icr 
error bars for integer q from —3 to 3. 



quantities A{q,<P,L) and F{q,<P,L) includes many more 
points for the scheme B as for the other two methods, 
which makes the fitting much more reliable for scheme B. 
For the other two schemes the actual number of consid- 
ered data points of course strongly depends on the number 
of factors of the system size N . Further, schemes A and 
B yield different values for the quantities A{q, <P, L) and 
F{q, (p, L) because scheme B considers all possible box ori- 
gins in contrast to scheme A. The biggest differences are 
observed for large box sizes L because in this case only 
very few box origins are taken into account in scheme A 
as visualized in Fig. [21 The reason is that the number of 
considered box origins has a crucial influence on the dis- 
tribution of the box probabilities /!;,(<?, L). This becomes 
clearly visible by plotting the distribution with and with- 
out an explicit average in Fig. [71 While the distributions 
for small box sizes L are quite similar, there are huge 
differences already for medium-sized boxes. Even for the 
small box sizes the distribution is not as smooth as the 
results averaged over all box origins. As we apply peri- 
odic boundaries there is no preferable choice for the box 
origins. Hence, the average over all box origins is a bet- 
ter estimate for the real box probability distribution. For 
the method C, which also includes an average over all box 
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Fig. 7: Distribution of box probabilities for the scheme A 
and scheme B for the wave function considered in Fig. |4] 
with linear box sizes L = 15 (top) and i = 30 (bottom). 



origins, the data points are identical to those of the parti- 
tioning scheme B but only box sizes for integer values of 
1/e are taken into account. Hence, the results for method 
B and method C are often rather close, but the accuracy 
of the fit for C is not as good as for B. 

Having a look at the singularity spectra /(ctq) we ob- 
serve that the results are rather close for all three schemes. 
However, the error bars are significantly larger for scheme 
A than for scheme B especially for the system size N — 
360. This is not surprising because in B many more data 
points are considered for the least squares fit and the av- 
erage over all box origins leads to a better estimate of the 
box probabilities fj,b- Comparing schemes B and C we ob- 
serve that the errors are rather similar, i.e., most of the 
improvement of the fits is due to the averaging over all 
different box origins. For the system with N = 350 the 
error bars in Fig. [6] are rather close for all three schemes 
especially for negative values of q and only a slight advan- 
tage of scheme B can be found. At first view this seems 
surprising because for this system even fewer box sizes 
are taken into account for scheme A and C. The reason 
is that by considering box sizes from L = 10 to N/2 in 
scheme B large boxes are strongly overrepresented in the 
least squares fits. This leads to a strong contribution of the 
non-linear behavior of A{q, L) and F{q, ^, L) for large 
L and negative q to the fit in B (cp. Fig. [5]) while for A 
and C almost no data points are considered in this range. 



5.2 Influence of range of box sizes L 

Further, we also investigated the influence of the consid- 
ered range of box sizes. As mentioned in Sec. [31 the ex- 
ponents should be computed from the scaling behavior 
of preferably small boxes and by considering the range 
1 <C i to properly measure the multifractal proper- 
ties. The first point is not properly considered yet because 
especially in scheme B larger box sizes L dominate the fit 
as mentioned above. In order to neglect the influence due 
to the different number of considered box sizes for the fits. 
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Fig. 8: Singularity spectrum f{aq) obtained using scheme 
A or scheme B with different ranges of box sizes L for the 
wave function in Fig. 2] with ~5 < g < 5 and la error bars 
for integer q from —4 to 4 (top). Corresponding squares 
of the correlation coefficients for the singularity strength 
aq and the fractal dimension f{aq) (center and bottom). 



we study the singularity spectrum for an identical number 
of box sizes for the different schemes. This means that we 
compute the singularity spectrum for the two-dimensional 
wave function in Fig. [T] by using exactly 15 different box 
sizes, i.e., for scheme A we have taken all box sizes with 
integer ratios N/L in the range from L = 10 to N/2 into 
account and for scheme B we used on the one hand L = 10 
to 24 and on the other hand L = 5 to 19. In the latter case 
we do not start with imin = 2 to assure i 1 and, there- 
fore, to omit the deviations from the power-law behavior 
which can be observed in Figs. [Hand [5] for small box sizes 
and negative values of q. This deviations are most likely 
caused by numerical inaccuracies which are enhanced for 
negative powers of /.tf, and by summing over many boxes. 
The results in Fig. [5] again show that the error bars are 
significantly smaller for scheme B. We can also observe 
clear differences in the singularity spectrum especially for 
negative values of q for all three considered setups. 
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We also plot in Fig. [5] the squares of the correlation 
coefficients obtained from the least squares fit, which is a 
measure for the amount of variability in the data that is 
accounted for by the fitting curve. Values of = 1 mean 
that the data can be fitted perfectly by a straight line, i.e., 
in this case we can identify a power-law behavior. While 
the numerical results in Fig.[8]show that for the singularity 
strength aq the quantity I for — 5 < g < 5, we 

observe smaller values of for large positive and negative 
q for the least squares fits of the fractal dimension fq in 
agreement with the increase of the fit errors for these q 
values. The results also depend on the used scheme and 
range of box sizes, but again scheme B yields the best 
results. Using scheme B with box sizes from L = 5 to 19 
the scaling behavior can be well described by a power-law 
for — 4 < g < 3 and 1 < < 2.8 respectively. 



5.3 Three-dimensional wave functions 

Finally, we have performed the same calculations also for 
the three-dimensional Anderson model for the disorder 
strength w — 16.5 at the mobility edge and a linear sys- 
tem size N = 60. At the mobility edge in three dimen- 
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sions the wave functions are known to be multifractals 
[3j. The results for the scaling of the quantities A{q,(!>,L) 
and F{q,(l>,L) are shown in Fig. Inland the correspond- 
ing singularity spectrum is shown in Fig. 1101 Again both 
quantities, A{q, <P, L) and F{q, <l>, L), show a power-law be- 
havior and yield qualitatively similar results as in two di- 
mensions for the partitioning schemes A and B. However, 
for the singularity spectrum we can observe a clear differ- 
ence. While the results for the schemes B and C are nearly 
identical, the singularity spectrum obtained by the parti- 
tioning scheme A shows significant differences for negative 
values of q (large aq). In Fig. [5] we see that the average 
values of A{q,<l>,L) and F{q,<P,L) are already quite dif- 
ferent for negative values of q. This is probably related to 
the rather small system size. Hence, the distribution of the 
box probabilities is already rather different for relatively 
small box sizes. Further, due to the identical distribution 
of box probabilities for the partitioning schemes B and 
C, the only differences are the different number of consid- 
ered data points for the least squares fit and hence the 
singularity spectra for both schemes are very similar. 

All results show that scheme B significantly reduces 
the fit errors. This can be expected to ultimately lead to 
a better convergence of the ensemble average. To demon- 
strate this, we also show the singularity spectrum for wave 
functions closest to E = for ten different disorder real- 
izations and the corresponding fit errors for the schemes A 
and B in Fig.[TT] For moderate values of q scheme B shows 
very small fluctuations of the scaling exponents, which are 
in all cases smaller than for scheme A. We also present the 
corresponding ensemble average from these ten wave func- 
tions in Fig. [T^l which is obtained from Jj 
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Fig. 9: Same as Fig. 01 but for a wave function of the three- 
dimensional Anderson Hamiltonian for the linear system 
size = 60 and the disorder strength w — 16.5. The least 
squares fits consider all box sizes in the range 5 < L < N/ 4 
marked by the two vertical dashed lines. 
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Fig. 10: Singularity spectrum f{aq) for the three different 
methods for the system in Fig. [9] with considered box sizes 
from L = 10 to N/2. The plot shows la error bars for 
integer q from —3 to 3. 
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Fig. 11: Smgularity spectrum f{aq) for wave functions at 
E = for 10 different disorder realisations of the three- 
dimensional Anderson Hamiltonian for the system size 
A'' = 60 and the disorder strength w — 16.5 computed 
with scheme B with 5 < L < 15 (top) and scheme A 
with 5 < L < 30 (bottom). The plot shows la error bars 
for integer q from —2 to 2, and identical disorder realiza- 
tions/wave functions are used for schemes A and B. 

The error bars indicate the standard deviation due to the 
least squares fit and are not due to the fluctuations of the 
points for different disorder realizations. Wc also show the 
squares r'j^^^ of the correlation coefficients of the fractal 
dimension for both methods. The standard deviations are 
quite small for both methods although we average only 
over ten wave functions. However, method B shows a sig- 
nificantly better performance for large positive and nega- 
tive q. 

We also like to point out that in literature different 
approaches for the determination of averages and errors 
in multifractal analysis can be found. On the one hand 
the exponents aq and fq were determined from a least 
squares fit according to Eq. pT|) (1|16) , and on the other 
hand they were computed from the average of aq{<l>) and 
/q(<?) obtained for the single wave functions |17I18| . The 
first method yields the ensemble average and the second 
method corresponds to the so-called typical average of the 
multifractal properties of the wave functions. For reasons 
of comparison we show in Fig.[l3]the results for the typical 
average of the ten disorder realizations in Fig. 1111 The 
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Fig. 12: Ensemble average of the singularity spectrum 
(solid lines) and the correlation coefficients ''^(q) of the 
fractal dimensions for the wave functions in Fig. [TTJ The 
error bars reflect only the accuracy of the least squares fit 
for q from —1.5 to 3. We check whether the results fulfill 
the symmetry relation in Eq. (jl2l) by plotting the quantity 
f{2d — aq) — d + aq (dashed lines). 




scheme B 



1 2 3 4 5 6 7 

average singularity strength (a^) 

Fig. 13: Typical average of the singularity spectra (solid 
lines) and the averaged correlation coefficients of 
the fractal dimensions for the wave functions in Fig. [TT] 
The error bars reflect the standard deviation from the av- 
eraged exponents (aq) and (fq) for q from —2 to 2. We 
again check whether the results fulfill the symmetry rela- 
tion in Eq. ([T^ (dashed lines). 



results for scheme A and B are quite close but the error 
bars for scheme B are again significantly smaller. The main 
differences between the ensemble and the typical average 
occur for large positive and negative q because for these 
cases the typical average is dominated by the behavior of 
a single wave function |14) . 
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Further, in literature the error bars denote often differ- 
ent errors of the method like e.g. the standard deviation 
of the least squares fit 5 19 or the standard deviation 
from averaging over the exponents obtained for different 
disorder realizations [17120] . More advanced analyses use 
error propagation and consider correlations for the data 
points in the least squares fits according to Eq. (ITlT) [TiT. 

In addition, we compare the resulting singularity spec- 
trum of the ensemble average in Fig. [T^] and of the typical 
average in Fig. [13] to the symmetry relation 

f{a,)^fi2d-a,)-d + a, (12) 

proposed by Mirlin et al. [H]. We observe that this sym- 
metry relation is fulfilled well for the ensemble and the 
typical average for positive q. However, for large Uq the 
ensemble averaged data show some differences from the 
symmetry relations, and we find considerable deviations 
for the typical average. Further, the proposed upper bound 
aq < 2d is satisfied neither for the ensemble average nor 
for the typical average [21, although the results of scheme 
B for the ensemble average in Fig. [12] are quite close to 
this upper bound within the error margins. These singu- 
larity strengths correspond to large negative values of q 
and, therefore, small errors are strongly enhanced in this 
region and the fit quality is decreasing (cp. also Fig. (S]). 
This has been found before in many other numerical stud- 
ies. For instance, the results by Rodriguez et al. show that 
the bound aq < 2d is approached for larger system sizes 
[16122] . However, the comparison between the numerical 
results for scheme A with this symmetry relation shows 
bigger differences. The observations for the symmetry re- 
lation and for the ensemble average in Fig. [12] support the 
finding that scheme B can help to improve the ensemble 
average in the multifractal analysis. 



6 Conclusion 

We have compared different partitioning schemes in the 
box-counting algorithm in order to improve the evaluation 
of Eqs. (O and ^ for a single wave function. We believe 
that the partitioning scheme considering unrestricted ra- 
tios 1/e of the system size and the box size as well as 
an average over all box origins is better suited than the 
best method proposed by Rodriguez et al. |T]. This con- 
clusion is based in particular on the observation that our 
method leads to smaller error bounds in the least squares 
fits and on the fact that it is not restricted to system sizes 
with preferable large numbers of factors. Our scheme is 
not meant to improve the ensemble average procedure and 
therefore it does not interfere with the established error 
estimates of this average. We improve the evaluation of 
Eqs. ^ and (jHl) by avoiding any bias in the choice of 
the origin of the boxes and by consequently averaging in 
Eq. (|10p over all possible choices. We found that the av- 
erage over different box origins yields a better estimate of 
the distribution for the box probability already for mod- 
erately large box sizes L. Due to the thus improved values 
of A{q,(l>,L) or F{q,'P,L), a better fit of the singularity 



strength aq and the fractal dimension fq becomes feasi- 
ble. Further, we showed that only small box sizes should 
be considered for the fit. Additionally, our scheme allows 
us to consider many more box sizes L, what also helps to 
improve the fits of aq and fq. 

Rodriguez et al. also studied the multifractal proper- 
ties of the Anderson Hamiltonian with system-size scal- 
ing although it is computationally more demanding [6114] . 
They found that system-size scaling yields better fit re- 
sults than the box-scaling approach. However, they com- 
pared their results with that of partitioning scheme A. 
Hence, a comparison with the better partitioning scheme 
B is necessary to validate the better performance of the 
system-size scaling approach. But we expect that also for 
the system-size scaling scheme B is superior in compar- 
ison with scheme A, not only because of the averaging 
over the box origins, which yields more accurate values of 
A{q,(l>,L) or F{q,<P,L), but also because there is no re- 
striction for the system sizes: in A again one needs a value 
of N which factorizes into many factors while in B every 
N can be included in the evaluation. 
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